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Abstract 

We analyze a class of exact type II solutions of the Robinson-Trautman family 
which contain pure radiation and (possibly) a cosmological constant. It is shown 
that these spacetimes exist for any sufficiently smooth initial data, and that they 
approach the spherically symmetric Vaidya- (anti-)de Sitter metric. We also inves- 
tigate extensions of the metric, and we demonstrate that their order of smoothness 
is in general only finite. Some applications of the results are outlined. 

PACS: 04.30.-w, 04.20.Jb, 04.20.Ex 

1 Introduction 

The classic Vaidya metric [1-4] (see also [5,6] followed by reprints of the original Vaidya 
papers) is a spherically symmetric type D solution of the Einstein equations in the presence 
of pure radiation matter field which propagates at the speed of light. In various contexts 
this "null dust" may be interpreted as high-frequency electromagnetic or gravitational 
waves, incoherent superposition of aligned waves with random phases and polarisations, 
or as massless scalar particles or neutrinos. The Vaidya solution depends on an arbitrary 
"mass function" m(u) of the retarded time u which characterises the profile of the pure 
radiation (it is a "retarded mass" measured at conformal infinity). Various sandwiches and 
shells of null matter can thus be constructed that are bounded either by flat (m = 0) or 
Schwarzschild-like (m = const 7^ 0) vacuum regions. Due to this property such solutions 
have been extensively used as models of spherically symmetric gravitational collapse of a 
star, as an exterior solution describing objects consisting of heat-conducting matter, as 
an interesting toy model for investigation of singularities and their possible removal by 
quantum effects, for studies of various formulations of the cosmic censorship conjecture 
on both classical and quantum level, process of black-hole evaporation, and for other 
purposes (see, e.g., [7-16] for more details and related references). 



In fact, the Vaidya spacetime belongs to a large Robinson-Trautman class of expand- 
ing nontwisting solutions [4,17,18]. Various aspects of this family have been studied in the 
last two decades. In particular, the existence, asymptotic behaviour and global structure 
of vacuum Robinson-Trautman spacetimes of type II with spherical topology were inves- 
tigated [19-28], most recently in the works of Chrusciel and Singleton [29-31]. In these 
rigorous studies, which were based on the analysis of solutions to the nonlinear Robinson- 
Trautman equation for generic, arbitrarily strong smooth initial data, the spacetimes were 
shown to exist globally for all positive retarded times, and to converge asymptotically to 
a corresponding Schwarzschild metric. Interestingly, extension across the "Schwarzschild- 
like" event horizon can only be made with a finite order of smoothness. Subsequently, 
these results were generalized in [32, 33] to the Robinson-Trautman vacuum spacetimes 
which admit a nonvanishing cosmological constant A. It was demonstrated that these 
cosmological solutions settle down exponentially fast to a Schwarzschild-(anti-)de Sitter 
solution at large times u. In certain cases the interior of a Schwarzschild-de Sitter black 
hole can be joined to an "external" cosmological Robinson-Trautman region across the 
horizon with a higher order of smoothness than in the corresponding case with A = 0. For 
the extreme value 9Am 2 = 1, the extension is smooth but not analytic (and not unique). 
The models with A > also exhibit explicitly the cosmic no-hair conjecture under the 
presence of gravitational waves. On the other hand, when A < the smoothness of such 
an extension is lower. 

Our aim here is to further extend the Chrusciel-Singleton analysis of the Robinson- 
Trautman vacuum equation by including matter, namely pure radiation. It was argued 
already by Bicak and Perjes [34] that with A = such spacetimes should generically 
approach the Vaidya metric asymptotically. We will analyze this problem in more detail, 
including also the possibility of A ^ in which case the Robinson-Trautman spacetimes 
containing pure radiation can be shown to approach the radiating Vaidya-(anti-)de Sitter 
metric. 



2 The metric and field equations 

In standard coordinates the Robinson-Trautman metric has the form [4,18,35] 

(m A \ r 2 — 

K - 2r (In P) u - 2 r 2 du 2 - 2dudr + 2— d(d( , (1) 
r 3 / P z 

where K = A(lnP) with A = 2P 2 d^d^ being the Gaussian curvature of the 2-surfaces 
2P~ 2 d(d(, m(u) is the mass function, and A is the cosmological constant. When the 
function P(u, (, () satisfies the fourth-order Robinson-Trautman field equation 

AK + 12m(lnP), u - 4m, u = 2nn 2 , (2) 

the metric describes a spacetime (generally of the Petrov type II) filled with pure radiation 
field Tju, = n 2 (u, (, () r~ 2 k^k V) where k = d r is aligned along the degenerate principal null 
direction (we use the convention G^ u + Ag^ u = kT^ u ). In particular, vacuum Robinson- 
Trautman spacetimes are given by n = 0, in which case m can be set to a constant by a 
suitable coordinate transformation [4]. Vacuum spacetimes (0) — possibly with a nonvan- 
ishing A — thus satisfy the equation 12 m (lnP) )U = —A A". These include the spherically 
symmetric Schwarzschild-(anti-)de Sitter solution which corresponds to P = 1 + |CC 
Indeed, replacing the complex stereographic coordinate ( with angular coordinates by 
( = v/2e^tan(fl/2), we obtain 2P ~ 2 dCdC = d6 2 + sin 2 fld0 2 , and K = A ln(P ) = 1. 
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Figure 1: Schematic conformal diagrams of the Robinson-Trautman exact spacetimes 
which exist for any smooth initial data prescribed on uq. Pure radiation field is present in 
the shaded region u < 0. Near u = the solutions approach the Vaidya metric, and can 
be extended to flat Minkowski region u > 0. Thick line indicates the curvature singularity 
at r = whereas double line represents future conformal infinity X + at r = oo (A = is 
assumed). The global structure depends on the value of the parameter fi of the linear mass 
function (jHJ): left diagram corresponds to fi > 1/16, the right one applies when /i < 1/16. 



Here we will restrict ourselves to nonvacuum cases for which the dependence of the 
mass function m(u) on the null coordinate u is only caused by a homogeneous pure 
radiation with the density n 2 (u)r~ 2 . When the mass function m(u) is decreasing, the 
field equation (j2J) can be naturally split into the following pair, 

AK + 12m(u) (lnP), u = 0, (3) 
—2m(u) tU = Kn 2 {u) . (4) 

In fact, it was demonstrated in [34] that such a separation can always be achieved using 
the coordinate freedom. It is then possible to reformulate equation (jHJ) by introducing 
a M-dependent family of smooth 2- metrics g a b on the submanifold r = const, u = const, 
such that g a i, = f(u, (, C)~ 2 9abi where g° b (C, C) is the metric on a 2-dimensional sphere S 2 . 
Since g a b is of the form 2P~ 2 d(d( in our case, we can write 



and equation (jHJ) becomes 



P = fP 0i P = 1 + |CC, (5) 
df 1 



du Ylmiu 



fAK, (6) 



where A is the Laplace operator associated with the metric g a b- Denoting A and K = 1 
as the corresponding quantities related to g® b , we obtain 



A = / 2 A , K = / 2 (1 + A (ln/)) 



(7) 



3 Linear mass function 



Let us first consider the simplest choice of m(u) which, in fact, has been widely used in 
literature (see e.g. [7,9,36]): we will assume that the mass function is a linearly decreasing 
positive function 

m{u) = /i = const > , (8) 

on the interval [uo, 0]. Notice that for (jHJ) the pure radiation field is uniform because equa- 
tion (jlj) implies n = k = const, independent of the retarded time u. The constant 
value «o < localises an initial null hypersurface (that extends between the curvature 
singularity at r = and the conformal infinity r = oo) on which an arbitrary sufficiently 
smooth initial data given by the function 

/o(C,C)=/(« = «o,C,C) , (9) 

are prescribed, see Fig. [T] 



3.1 Existence of the solutions 

Now, the idea is to employ the Chrusciel-Singleton results [29-31] concerning the analysis 
of the Robinson-Trautman vacuum equation, in particular the existence and asymptotic 
behaviour of its solutions. In the vacuum case m in equation (JHJ) is constant, and the 
solution f(u,(,C) °f t ne characteristic initial value problem Q exists and is unique (in 
spite of the singularity at r = 0). In the presence of pure radiation given by (jHJ) it is 
possible to "eliminate" the variable mass function from the Robinson-Trautman field 
equation © mathematically by a simple reparametrisation 

u = -/i" 1 ln(-w) , (10) 

cf. [34]. Indeed, equation (JHJ) is then converted to 

df 1 

where /(«, C, C) = /(«(«), CO, K = f 2 (l + A ln(/)), and A = f 2 A . Notice that the 
transformation (fTU|) moves the hypersurface u = 0, on which the mass function m(u) 
reaches zero, to u = +oo. 

Chrusciel [30] derived the following asymptotic expansion (as u — > oo) for the func- 
tion / satisfying the evolution equation (fTTj) for any smooth initial data f = f on 
u = — iT 1 ln(— u ), namely 

f = l + /i,oe- a + / 2 ,oe- 4fi + -- + /i4,oe-* 

+/i5,iie- 30 " + / 15 ,oe- 30 " + --- (12) 

oo Ni 
i=0 j=0 

where f\j are smooth functions on S 2 such that fij = for j > 0, i < 14. The function / 
thus exists and converges exponentially fast to 1, which means physically that the radia- 
tive Robinson-Trautman vacuum spacetimes approach asymptotically the Schwarzschild- 
(anti-)de Sitter solution as u — > oo, see relation (j3J). In our case of pure radiation field (JBJ) 



we employ the transformation (fTUj) on expression (jl2j) to obtain the following asymptotic 
expansion of / as u — > 0_, 

/ = 1 + / li0 {-uf/» + / 2)0 (- U ) 4/ " + • • • + /l4,0 (-«) 28/M 

-Ai-Vis,! ln(-u) (-«) 30/ ^ + /i 5 ,o (-«) 30//t + ■ • • (13) 

oo Ni 

= EE/«[-^" l M-«)l i (-«) ai/M - 

As a result, for the initial data Q the Robinson-Trautman type II spacetimes which 
contain uniform pure radiation field with the linear mass function (JHJ do exist in the 
whole region u < u < 0. It is also obvious that the function / approaches 1 as u —* 0_ 
(where also m(u) — > 0) according to ([13)1 . In other words, these spacetimes approach the 
spherically symmetric Vaidya-(anti-)de Sitter metric near w = 0. 

The global structure of such spacetimes is schematically indicated on Fig. ^ In fact, 
there are two possibly different conformal diagrams depending on the value of fx: for 
fx > 1/16 there is a white hole singularity at r = 0, for fx < 1/16 there is also a naked 
singularity, see e.g. [9,13,16,36] for more details. At u = all of the mass m(u) is radiated 
away, and we can attach Minkowski space (de Sitter space when A > 0, anti-de Sitter when 
A < 0; the presence of the cosmological constant would change the character of conformal 
infinity I which would become spacelike or timelike, respectively) in the region u > along 
the hypersurface u — 0. We will now investigate the smoothness of such an extension. 

3.2 Extension of the metric across u = 

It follows from (|13)1 that the smoothness of / on u = is only finite. Depending on the 
value of fx two different cases have to be discussed separately: 2/ fx is an integer, and 2/ \x 
is a real non-integer positive number. 

When 2 l\i is an integer then due to the presence of the ln(— u) term associated with 
/i5,i 7^ the function / is of the class C^ 30 ^ -1 . For ft very small, the integer number 
(30///) — 1 is large so that / becomes smoothly extendable to 1 across u = as fi ~^ 0. 
This represents a naked-singularity Robinson-Trautman spacetime (see the right part of 
Fig. HJ unless fi = which gives flat space everywhere. In the limiting case \l — 1/16 the 
function / is of the class C 479 . For the (white hole) Robinson-Trautman spacetimes given 
by /i > 1/16 the smoothness is lower. However, it is always at least C 14 because fi < 2 in 
this case. 

In the generic case when 2 /fx is not an integer the function / is only of the class C^ 2 ^\ 
where the symbol {x} denotes the largest integer smaller than x. Again, with fi — > the 
function / becomes smoothly extendable. For fx < 1/16 the function / is at least of the 
class C 32 , for fx > 2 it is not even C l but it remains continuous. 

To investigate further the smoothness of the metric when approaching the hypersur- 
face u = 0_ which is the analogue of the Schmidt-Tod boundary of vacuum Robinson- 
Trautman spacetimes [25, 30] we should consider the conformal picture using suitable 
double-null coordinates. Such Kruskal-type coordinates for the Vaidya solution with lin- 
ear mass function (jHJ) were introduced by Hiscock [7-9], see also [16,36], and we will use 
this transformation only to replace the coordinate r since the null coordinate u is already 
appropriate. Introducing a new coordinate w by 



we put the Robinson-Trautman metric with linear mass function into the form 




y 2 



2r + u + 2/2— dudw + 2— dCdC , (15) 



r 



p2 



where r(u,w). For the pure Vaidya metric characterized by / = 1 and Kq = 1 the first 
term vanishes identically so that the coordinates of (fTo"]) are indeed the Kruskal-type 
coordinates for the Vaidya spacetime with a linear mass function. 

The smoothness of a general Robinson-Trautman metric (JT5j) depends only on the 
smoothness of the metric coefficients g uu and (containing the function /) since the 
coefficient g uw tends to — r as u — > 0. The smoothness of g^ (for any finite r) and of K is 
the same as of /, see (J7|). The function f )U /f is evidently one order less smooth than /. 
Consequently, for 2/ ' \x being integer or non-integer number, the metric (JT5j) is of the class 
c Wri-2 or C {2ffi-x^ reS pectively We again observe that the spacetimes approaching 
the linear Vaidya metric with naked singularity (i.e., for small values of the parameter /i) 
possess higher order of smoothness at u = 0. 

One might be worried about the invariance of our results, namely with respect to a 
rescaling of the null coordinate u(u) leading to a different smoothness of the function / 
and of the metric. In order to change the smoothness on the hypersurface u = such 
rescaling must have a singular character there. But this would lead to a degeneracy of the 
metric coefficient g^ w of the Vaidya metric, which is forbidden. Consequently, the above 
results are in this sense unique. 

We would like to obtain analogous results concerning smoothness of the extension also 
for a non-zero value of the cosmological constant A. Unfortunately, as far as we know, 
there is no explicit transformation of the Vaidya-de Sitter metric to the Kruskal-type 
coordinates even for the linear mass function (contrary to the Schwarzschild-de Sitter 
case [33]). However, it is possible to start with the Vaidya-de Sitter metric 

ds 2 = -h(u,r)du 2 -2dwdr + r 2 dVt 2 , (16) 

where h(u,r) = 1 + 2/iwr~ 1 — |r 2 , and perform a coordinate transformation 

dw = gdu + 2^dr , (17) 
h 

where g(u,r) is some function. We arrive at the double-null form for the metric 

ds 2 = --dudw + r 2 (u,w)dtt 2 . (18) 
9 

Of course, we have to ensure that dw in (|T7jl is a differential of the coordinate w. The 
integrability condition (d 2 w = 0) gives the following quasilinear PDE, 

h — - 2h— + A-g = , (19) 
or ou r 

for the undetermined function g, which is difficult to solve analytically. The method of 
characteristic curves leads to the first-order ODE of the Abel type which has not yet 
been solved, but the existence of its solution is guaranteed. [It is possible to apply the 
perturbative approach starting from the solvable case of the de Sitter metric (/i = 0) and 



then linearise the PDE in the parameter /x. The result, however, can not be presented 
in a useful closed form.] For our purposes it suffices to use a general argumentation: 
the coordinate u is already suitably compactified and we are only determining the com- 
plementary null coordinate w to obtain the Vaidya-de Sitter metric in the Kruskal-type 
coordinates (which is smooth on u = 0). The corresponding Robinson-Trautman metric 
in these coordinates differs only by the term g uu {u, r, (, Q du 2 (which is absent in the 
Vaidya-de Sitter case in the double null coordinates), and by a different metric coefficient 
(jf££ = r 2 / _2 P ~ 2 , where r(u,w) is finite and smooth when approaching the hypersurface 
u = 0. The smoothness is thus not affected by the specific transformation (|17j) and it 
is the same as for the vanishing cosmological constant. This is different from vacuum 
spacetimes with m = const 7^ studied in [32,33] because in the present case m — > near 
u — 0, and the influence of A on the smoothness becomes negligible. 



4 General mass function 

The results obtained above can be considerably generalized. Inspired by a similar idea 
outlined in [34] we may consider a reparametrisation on the null coordinate u by 

u = 7 (u) , (20) 

where 7 is an arbitrary continuous strictly monotonous function. We start with the 
evolution equation (|TT| for which the existence and uniqueness of solutions has been 
proven, and their general asymptotic behaviour ()12j) has been demonstrated. Now, by 
applying the substitution (J2H|) in equation (JTTJ) we obtain 

(where the dot denotes a differentiation) which is the evolution equation for the function 
f(u, (,(,)■ This is exactly the Robinson-Trautman equation (JHJ) for the mass function 

m{u) = -L . (22) 

For a given smooth initial data on w there thus exists the Robinson-Trautman spacetime 
([TJ, including the cosmological constant A, with the mass function (J22j) . To obtain a 
positive mass we consider a growing function 7(1*). Considering (jlj) this corresponds to a 
universe filled with homogeneous pure radiation 



2 7 



n\u) = -^- 2 . (23) 
k 7^ 

For consistency the function 7 must be convex. An asymptotic behaviour of the function 
/ as j(u) — > 00 is easily obtained from the expansion (|T2|) by substituting relation (f2~Uj) . 

In particular, the linear mass function (jHJ) discussed above is a special case of (J2*2~j) for 
the transformation (j20j) of the form (JTUJl . More general explicit solutions can be obtained, 
e.g., by considering the power function 

7 (u) = (-u)~ p , p > , (24) 

which gives 

m (u) = - (-u) 1+p , n 2 (u) = 2ijp+V) f- u )P . (25) 

p Kp 



Both functions m and n approach zero as u — > 0. Due to the theorems mentioned above, 
there exist Robinson-Trautman type II spacetimes in the region u < which approach 
the spherically symmetric Vaidya-(anti-)de Sitter metric as u — > 0_ with the mass func- 
tion and pure radiation given by (J25J) . The asymptotic behaviour of such solutions is 
determined by expression © with 

oo Ni 

f = i + E E /y (-«)"* ex P [-^(-^)~ p ] - (26) 

i=l j=0 

where /jj = for j > if % < 14. Interestingly, the function / is now smooth on u = 
for any power coefficient p, but this still does not guarantee that the extension into flat 
region u > is analytic (see [33] for a similar situation). 
Another simple explicit choice is 

7(u) = —M~ l In [sinh(-w)] , M > , (27) 

which implies (see also [16]) 

2M 

m(u) = M tanh(-w) , n 2 (u) = -g— . (28) 

k cosh M 

In the region u < the mass function monotonically decreases from M to zero, while the 
pure radiation field grows from zero to the value 2M/ k as u —>■ 0. Let us note that in this 
case the integrated radiation density is finite on the interval (— oo, 0), f®^ n 2 (u) = 2M/k. 
The expansion near u = 0_ is 

oo Ni 

f = 1 + E E hi (-M- 1 In [smh(-u)]y sinh 2l / M (- M ) . (29) 

i=l 3=0 

If 2/M is an integer then the function / belongs to the class C( 3 °/ M ) _1 ; otherwise it is of 
the class C {2/M} . 



5 Possible modifications and applications 

The Robinson-Trautman pure radiation solutions in the region uq < u < approach- 
ing the Vaidya metric near u — 0, which can be extended (albeit non-smoothly) to flat 
Minkowski space in the region u > as in Fig. may be used for construction of various 
models of radiative spacetimes. For example, it is natural to further extend the solution 
"backwards" into the region u\ < u < uq by the Robinson-Trautman vacuum solution 
with a constant mass m = m(uo), such that the function / is continuous on u . This is 
shown in Fig. |21 In such a case the spacetime may describe the process of "evaporation" 
of a white hole (with a different character of the singularity at r = when fi < 1/16) with 
its mass decreasing from the value m to zero. Let us emphasize that the region u < w 
does not represent the Schwarzschild solution because the spacetime is not spherically 
symmetric there (/ 7^ 1). In fact, this is the region where the original Chrusciel theorems 
on the behaviour of the Robinson-Trautman vacuum spacetimes with constant mass ap- 
ply (cf. (JUJ), (JEJ))- However, the spacetime in this region can not be extended up to the 
past conformal infinity I~ because the metric function / diverges as u —>■ — 00. 

In the presence of the cosmological constant A one obtains a family of exact spacetimes 
that describe evaporation of a white hole in the (anti-)de Sitter universe. In this case the 



r = 
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Figure 2: Possible extensions of the Robinson-Trautman radiative spacetimes into the 
region u < uq. Pure radiation is present only in the shaded region, everywhere else it is a 
vacuum solution. For u G (ui,u ) the mass function is constant, m(u ) = —fiu , but the 
spacetime is not spherically symmetric — it is not the Schwarzschild solution (/i > 1/16 
on the left, /J, < 1/16 on the right). 



schematic conformal diagram on Fig.|2]has to be modified in such a way that for all values 
of u the conformal infinity 1 + becomes timelike (for A > 0) or spacelike (for A < 0). 

Another possible modification is to consider the "advanced" form of the spacetimes 
(which describes an ingoing flow) rather than the "retarded" form (corresponding to 
outgoing flow) employed above (see, e.g., [13] for more details). This time- reversed form 
is obtained formally by a simple substitution u — > — v in the metrics and corresponding 
functions. The Robinson-Trautman metric thus reads 

(777 A \ T 2 — 

K + 2r(ln P) v - 2 r 2 dv 2 + 2dvdr + 2— d(d( , (30) 
r 3 / P 2 

where m(v) is an increasing mass function in v G [0, vq] . This is joined with flat Minkowskian 
region v < 0, and extended to the region v > v o by the corresponding Robinson-Trautman- 
(anti-)de Sitter black hole vacuum solution, see Fig. El It is a non-spherical generalization 
of the gravitational collapse of a shell of null dust forming a naked singularity [9,37,38] - 
in these works the mass function was taken to be m(v) = fiv (with m(v) = for v < 0, 
and m(v) = M = fivo for v > vq). The metric function P is now given by P = fPo where 
/ is analogous to (fT3"j) . 

oo Ni 

/ = EE/« (-/x-Mn^V^, (31) 

i=0 j=0 

so that the smoothness of the metric on the boundary v — depends on the parameter \x. 
For v G (vq, Vi) the spacetime is vacuum but not spherically symmetric. The metric 
diverges as v — > oo. Our results can thus be interpreted in such a way that — at least 
within the Robinson-Trautman family of solutions — the model [9] of collapse to a naked 
shell-focusing singularity which is based on the spherically symmetric Vaidya metric is 
not stable against perturbations. 



Figure 3: Time- reversed version of Fig. ^represents the "advanced" form of the Robinson- 
Trautman spacetimes (|3~U|l which describes an ingoing flow of radiation. 

6 Concluding remarks 

In our contribution we have analyzed exact solutions of the Robinson-Trautman class 
which contain homogeneous pure radiation and a cosmo logical constant. This is a natural 
extension of previous works [19-34] on properties of vacuum spacetimes of this family. 
We have demonstrated that these solutions exist for any smooth initial data, and that 
they approach the spherically symmetric Vaidya-(anti-)de Sitter metric. It generalizes 
previous results according to which vacuum Robinson-Trautman spacetimes approach 
asymptotically the spherically symmetric Schwarzschild-(anti-)de Sitter metric. We have 
investigated extensions of these solutions into Minkowski region, and we have shown that 
its order of smoothness is in general only finite. Finally, we suggested some applications 
of the results. For example, it follows that the model of gravitational collapse of a shell 
of null dust diverges as v — > oo which indicates that investigations of such process based 
on the spherically symmetric Vaidya metric are, in fact, not stable against "non-linear 
perturbations" , at least within the Robinson-Trautman family of exact solutions. 
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